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Abstract. We present the results of a simulation of a wind-driven non-linear gravity wave breaking 

qq ■ on the surface of a white dwarf. The "wind" consists of H/He from an accreted envelope, and the 

t^j- \ simulation demonstrates that this breaking wave mechanism can produce a well-mixed layer of 

(SJ ■ H/He with C/O from the white dwarf above the surface. Material from this mixed layer may then be 

f*"*. ' transported throughout the accreted envelope by convection, which would enrich the C/O abundance 

' of the envelope as is expected from observations of novae. 
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INTRODUCTION 



2 ■ Classical novae result from the ignition (and subsequent explosive thermonuclear burn- 

ing) of a (~ 10 m) layer of hydrogen-rich material that has accreted from a main se- 
quence companion onto the surface of a white dwarf [jl|, |3], Observed abundances 
and explosion energies estimated from observations indicate that there must be signif- 
icant mixing of the heavier material of the C/O or O/Ne white dwarf into the lighter 
accreted material (H/He). This mixing is critical because otherwise hydrogen burning 
would be too slow to reproduce observed nova characteristics in outburst. Further, with- 
out this mixing it is difficult to understand the observed abundances of intermediate-mass 
nuclei in the ejecta. Accordingly, nova models must incorporate a mechanism that will 
dredge up the heavier white dwarf material [|5], and references therein] . 

A recently proposed mixing mechanism is the breaking of non-linear resonant gravity 
waves at the C/O surface [^j, [7], |5|]. The gravity waves, driven by the "wind" of accreted 
material, can break, forming a layer of well-mixed material. This mixed layer may then 
be transported upward by convection, thereby enriching the accreted material. Because 
the length scale of this mixed layer may be very small (much smaller than the length 
scale of convection), previous precursor simulations have not captured this effect. 

In this manuscript, we present a simulation of a wind-driven non-linear gravity wave 
breaking on the surface of a white dwarf. The simulation was performed with FLASH, 
a parallel, adaptive-mesh simulation code for the compressible, reactive flows found in 
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FIGURE 1. Results from the simulation of a wind-driven gravity wave. The left panel shows the 
potential energy of the wave vs. time. Also shown is the potential energy predicted by the linear theory. 
The right panel shows the mixed C/O mass per unit area vs. time. 



many astrophysical environments [|| ^p. This simulation is part of an ongoing study 
of this mechanism to assess its efficacy for mixing white dwarf material with envelope 
material. 



THEORY AND SIMULATION DETAILS 

This work applies the theory of gravity wave generation (originally developed for the 
air over water interface) [[H], and references therein] to a white dwarf with an accreted 
envelope. The control parameters for this model are G = g8 /U^ ax , which is the ratio of 
potential to kinetic energy of the wind (g is the constant gravitational acceleration, t/ max 
is the maximum wind speed, and 8 is the characteristic length scale of the wind profile) 
and the ratio of densities at the interface, r = p l /p 2 - For this simulation, the domain 
was 1.0 x 10 6 cm by 1.0 x 10 6 cm, g = 4.5 x 10 9 cm/s, and r = 10, with a density of 
the white dwarf material at the interface of 10 4 gm/cm 3 . The white dwarf and accreted 
envelope materials are modeled as simple y = 5/3 gases with the white dwarf material 
composed (by mass) of a 50/50 C/O mix and the accreted material composed of a 75/25 
H/He mix. 

The density and pressure profiles were obtained by integrating the equation of hydro- 
static equilibrium 

dy 

which for the case of a compressible, gamma-law gas gives 

i 

p = pJ 1 -(r-i)yf]" 1 (2) 



FIGURE 2. Images of the log of density with velocity vectors at earlier times during the simulation. 
The units of density are g/cm 3 and the length of the velocity arrows is proportional to the magnitude of 
the velocity, with a maximum of 2 x 10 8 cm/s. The left panel shows the initial conditions, and the right 
panel shows the simulation at t = 0.015 s. 
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Here P is the pressure at the interface and p i is the density immediately above or 
below the interface. The wave was created by forcing the interface to be sinusoidal and 
perturbing the pressure and adding a velocity via a prescription similar to linear theory. 
The wind profile is given by 



tf(y) = tw(i-«" 



-y/s 



(4) 



with 5 = 1.0 x 10 5 cm and U max = 2 x 10 8 cm/s. The boundary conditions were 
isothermal hydrostatic on the upper and lower boundaries and periodic on the sides. 
The simulation was performed with seven levels of adaptive mesh refinement, for an 
effective resolution of 512x512 zones. 



RESULTS 

Figure [j] shows the potential energy of the wave and the amount of mixed material during 
the course of the simulation. The left panel is a plot of potential energy vs. time showing 
the result from the simulation and the potential energy given by the linear theory as 
4.0 x 10 26 exp(140f). The growth rate in this expression (140) comes from the control 
parameter G, and 4.0 x 10 26 (erg) is the initial potential energy of the wave [|7p. The two 
potential energies agree reasonably well until the wave breaks at about 0.04 s. The right 
panel is a plot of mixed C/O mass per unit area vs. time. The figure demonstrates the 
dramatic increase in the amount of mixed material that occurs above t = 0.25, the point 
during the course of the simulation at which the wave begins to break. 



FIGURE 3. Images of the log of density with velocity vectors at later times during the simulation. The 
units of density are g/cm 3 and the length of the velocity arrows is proportional to the magnitude of the 
velocity, with a maximum of 2 x 10 8 cm/s. The left panel shows the simulation at t = 0.030 s, and the right 
panel shows the simulation at t = 0.045 s. 



Figure |2| shows gray scale images of the log of density with velocity vectors at early 
times in the simulation. The longest velocity arrows, at the top of the images, correspond 
to 2.0 x 10 8 cm/s, the maximum wind speed. The left panel shows the initial conditions, 
and the right panel shows the configuration at t = 0.015 s. Visible in the right panel is 
the development of Kelvin-Helmholtz instability on the upwind (left) side of the wave. 
Figure [| shows similar gray scale images of the log of density with velocity vectors at 
later times in the simulation. By t = 0.030 s (left panel), the wave has begun to break, 
and the images show vorticity downwind from the crest of the wave. At t = 0.045 s (right 
panel), there is substantial mixing and obvious vorticity. 



CONCLUSIONS 

The simulation presented in this manuscript demonstrates the proposed breaking wind- 
driven non-linear gravity wave mixing mechanism. The results show the development 
of well-mixed zone just above the surface of the white dwarf. This simulation will be 
one part of a study of this mixing mechanism investigating effects of wind profiles and 



speeds. Complete details of the study will appear in Alexakis, et al. QTTJ] . 

The expectation is that the breaking of non-linear gravity waves on the surface of the 
white dwarf will lead to a thin well-mixed layer of material that may then be transported 
throughout the the envelope by convection. This study investigating the mixing mecha- 
nism should provide quantitative information about the mixing rate that will allow for 
the development of subgrid models that may be applied to multidimensional convection 
simulations to study the enrichment of the envelope. A preliminary simulation of this 



kind is also presented in this volume yl2Q. 
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